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Constraints on the anisotropy of dark energy 
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If the equation of state of dark energy is anisotropic there will be additional quadrupole anisotropy 
in the cosmic microwave background induced by the time dependent anisotropic stress quantified 
in terms of Aw. Assuming that the entire amplitude of the observed quadrupole is due to this 
anisotropy, we conservatively impose a limit of |Atu| < 2.1 x 10~ 4 for any value of w > —1 assuming 
that fi m < 0.5. This is considerably tighter than that which comes from SNe. Stronger limits, upto 
a factor of 10, are possible for specific values of f2 m and w. Since we assume this component is 
uncorrelated with the stochastic component from inflation, we find that both the expectation value 
and the sample variance are increased. There no improvement in the likelihood of an anomalously 
low quadrupole as suggested by previous work on an elliptical universe. 

PACS numbers: 98.80.Cq, 98.80.Jk 



Introduction: There is now very strong evidence that 
expansion of the Universe is accelerating based on a 
combination of measurements of type la supernovae 
(SNe) [H, [1], anisotropies of the cosmic microwave back- 
ground (CMB) [1,0] and large scale structure Many 
different explanations for this phenomenon have been 
suggested with most of them requiring the existence of 
new physics at a scale comparable to the present hori- 
zon size ~ c/Hq (for example, refs. [1, 0])- One of these 
ideas is that there is some new energy component, known 
as dark energy, with an equation of state P = wp and 
w < —1/3 in order to achieve acceleration. 

Measurements of CMB anisotropies appear to show 
that there a number of anomalies in their detailed 
statistical properties [1, H QjJ (HI- These suggest 
that the anisotropics may not be compatible with an 
isotropic, Gaussian random field on the very largest 
scales. Given that the Fourier modes corresponding to 
these anisotropies crossed the horizon in the last few 
Hubble times, that is, when the dark energy has come 
to dominate the expansion of the Universe, it is tempt- 
ing to connect the two. 

The equation of state parameter w is often considered 
to be a function of time, but the overall pressure tensor, 
Pi° , is usually assumed to be isotropic, that is Pi 3 = 
wpdcbi where pdc is the density of the dark energy. In 
this letter we will consider the possibility that the dark 
energy is parameterized by 



Pi- 



Pie [w5 t 3 + Aw. t 



(1) 



where we impose the traceless condition Aw, 1 = 0. We 
will show that an anisotropic equation of state leads to an 
additional component in the quadrupole anisotropy and 
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hence the observed amplitude of the CMB quadrupole 
leads to a strong constraint on Awy which is much 
stronger than those which have been deduced from SNe 
data |ll [HE!. 

The st udy of anisotropic universes has a long his- 
tory [l5|, 0, E3i but in the past the possibility that 
the Universe is rotating, being modelled by a non-trivial 
Bianchi universe, was the main focus. In our work we 
will concentrate on models which only have time depen- 
dent and spatially homogeneous anisotropic stress and 
no vorticity. Moreover, it is only the dark energy com- 
ponent which is anisotropic allowing the evolution of the 
Universe to proceed in the standard way until the dark 
energy begins to dominate. 

An anisotropic equation of state is natural in the Elas- 
tic Dark Energy (EDE) model [H, [Hj] , whereby the dark 
energy component is comprised of some continuum fluid 
with a non-zero shear modulus, similar to the idea of 
a crystalline material. One particular manifestation of 
this idea could be a soap film type structure formed 
from domain walls at a low energy phase transition with 
w = —2/3 [2Cjj . However, the limits which we will derive 
apply to any particular manifestation of anisotropic dark 
energy, for example, as discussed in refs. [H (HI 

The reason for this is that, to first order, the 
anisotropic evolution due to Awi 3 and the evolution of 
initial metric perturbations decouple and can be com- 
puted independently, then added at the end. If we use 
the time dependent s pat ial metric 7^ (77) to describe the 
anisotropic evolution |2l| and hij(rj, x) to represent effect 
the metric perturbations created during inflation, then 
our approach works when py^l 3> \hij\. It should also 
yield qualitative information when 7^ ~ hij . 

Anisotropic universes: We will consider universes with 
metric 



ds z 



a 2 {-dry 2 



ji 3 (r])dx l dx 3 ] 



(2) 



for which V' 
1 d 



HSi 3 +0-^ and T° 



HSij + a a , where 



<7jj = 2 dtjlij an d "H = a' /a, with the derivative with re- 



2 



spect to conformal time, r\. The stress-energy content will 
be made up of dark matter (m) , which will be isotropic 
and prcssureless, and dark energy (de), where the pres- 
sure tensor of the dark energy is given by |T]). The Ein- 
stein and conservation equations can be written as [2l| 



m 1 



Pdc 



= 8nGa p t0 
= -3W(H 



W)pde 



■j l Awi 3 pde , 



8irGa 2 Awi J pd e 



(3) 

(4) 
(5) 



where p tot = p ln + p dc , & 2 = o-'i 3 o~j l and we have used 
a flat geometry. One might think that there would be a 
velocity generated , but this would only be created by the 
gradient of the anisotropic stress which is absent in our 
model. Hence, we can consistently set this to zero and 
ignore it. Assuming there is no primordial anisotropic 
stress, that is, initially = 0, these equations can be 
solved to give pdc and as functions of a if Aw^ is 
considered small. To order (Aw) 3 , we find that 

Pde(t ) 



Ale (a) 

Qj 3 (a) 
H 

where E(a) 



3(l+to) « 

3tt dc Awi j 



1 



3fi de Aw l i Aw;/G(a)) 



(O m /a 3 
F(a) = 

G{a) = 



-F(a) , 
Ode/a 3(1+w) ) 1/2 , 



db 



, b^ w E{b) ' 
a dbF(b) 



(6) 
(7) 

(8) 
(9) 



and fix is the density of species X relative to the present 
day critical density, with f2 m + Q^c = 1- The function 
F(a)/a 2 is the growth function for anisotropic stress. 

The anisotropic stress grows quickly during matter 
domination with F(a) « a( 3 / 2 )- 3, "/(0^ 2 ((3/2) - 3w)) 
and the rate of increase declines during the dark en- 
ergy dominated era. We have plotted d(logF)/(3/2 — 
w)d(loga) against a for a range of values for w — — 1, 
-2/3 and -1/3 for Q m = 0.3 in Fig.Q] The slopes start 
at 4.5, 3.5 and 2.5, respectively, and decrease with time; 
this happens later in models with more negative values of 
w, since dark energy domination happens later for fixed 

ri m . 

Calculation of temperature anisotropics: The CMB 
anisotropics due to the anisotropic equation of state are 
given by 

AT f Vo ■ ■ 

— (n) = -J (Tijfffv 1 d<q , (10) 

where t]q and rj vec are the conformal time at the present 
and recombination eras, respectively. Using the expres- 
sion above with <xy = "fikO~j k and converting the integra- 
tion variable to a (ignoring terms higher order in Aw), 
one finds that 



AT 



(n) = —Awijffh 3 J(Q, m ,w) . 



(11) 




FIG. 1: The power law slope of the function F(a) for w — — 1 
(solid line), —2/3 (dotted line), —1/3 (dashed line). In each 
case we have used Q m — 0.3. 



where 



j(n mj «;) = 3(i-n m ) 



da 



db 



a 4 E(a) J b 1+3w E(b) 



(12) 

and <x roc ~ 1/1090 is the scale factor at recombination. 
Calculation of J(f2 m , w) will be insensitive to the precise 
value of a rec used because it is dominated by the late 
time behaviour. We have plotted J(Sl m , w) against w for 
a range of values of f2 m in Fig. [21 It is clearly 0(1) for 
the range of values relevant to observations. 

It may appear that the anisotropic equation of state 
parameter, Awij, has 5 degrees of freedom since it is a 
traceless, symmetric, rank 3 matrix. However, one can 
rotate the coordinate system of the observation, such that 
n can be replaced by Rh. This corresponds to a diagonal- 
ization of Awij and hence it is sufficient to parameterize 
Awij = diag(Awi, Aw2, — (Aw± + Au^)) and the rota- 
tion R which defines the direction of the anisotropy on 
the celestial sphere. The limits we will compute will not 
depend on R. Hence, we can write 



AT 

-^-(ft) = -J(Sl m , w) ( Awi sin 2 9 cos 2 



+Au>2 sin 2 9 sin 2 — (Awi + AW2) cos 2 



(13) 



It is clear that the effect of the anisotropic equation 
of state is to modify the temperature quadrupole so that 



FIG. 2: The function J(f2 m ,u>) against w for fi m = 0.1 
(dashed line), 0.2 (dotted line), 0.3 (solid line), 0.4 (short- 
dashed line) and 0.5 (dotted-short dashed line). 



FIG. 3: The 2a upper limit on the anisotropy, Awn m , against 
w for Jl m = 0.1 (dashed line), 0.2 (dotted line), 0.3 (solid line), 
0.4 (short-dashed line) and 0.5 (dotted-short dashed line). 



the multipole coefficients are given by 



02,2 
012,1 

fl2,0 
0-2,-1 

0-2,-2 



i 2,0 
4,-1 ) 

l 2,-2 ~ 



— J(Q ia ,w)(Aw 1 - Aw 2 ) ■ 
15 



J(Ct m ,w)(A Wl + Aw 2 ) 



-J(tt m ,w)(Aw 1 -Aw 2 ), (14) 




where a 2m is the stochastic quadrupole generated by ini- 
tial metric perturbations, which need not have a diagonal 
covariance matrix (for example, see ref. [24[ ). Using this, 
and assuming the two components arc uncorrclated, we 
can compute C 2 , the quadrupole coefficient of the angular 
power spectrum, whose mean value will be C 2 = C^ + C^. 
The deterministic power due to the anisotropic dark en- 
ergy is 



C£ = ^{J(n m ,w)] 2 (Aw) 2 



75 

with (Aw) 2 = Aw i j Awj i = 2(Awj 



Aw 2 



(15) 



Note that the effect on the average power spectrum is ad- 
ditive. Rather than having the usual \ 2 distribution with 
5 degrees of freedom, the quadrupole likelihood now has 
a non-central % 2 distribution. Writing the ratio of power 
between the deterministic and stochastic components as 



C 2 /C\, the variance of C 2 is given by 
(AC 2 f = l%±^ Cj. 



5(1 + aY 



(16) 



In Fig. 0] we illustrate the quadrupole likelihood func- 
tion for various non-zero a. In summary, the expecta- 
tion value of C 2 will increase and the fractional cosmic 
variance will decrease for a non-zero deterministic com- 
ponent. 

Limits on the anisotropy of dark energy: One can 
obtain a conservative constraint on Aw by assuming 
that the entire quadrupole amplitude comes from the 
anisotropic effect. The result from WMAP5 is that 
C 2 < 4.0 x lfT 10 at 95% confidence [H, H|, which im- 
plies a constraint of \Aw\ < 1.1 x 10~ 4 for w = —1, 
6.0 x 10- 5 for w = -2/3 and 2.2 x 10~ 5 for w = -1/3, all 
for Q m = 0.3. The upper limit, Au>n m = |Aiu|, for other 
values of w and f2 m is presented in Fig. [3] Making the 
conservative assumption that f2 m < 0.5 we deduce that 
| Aw | < 2.1 x 10 -4 . These limits are significantly stronger 
than those which come from type la SNe [H, Q which 
are typically |Ato| ~ 0.1. 

Conclusions: We have found that the consequences of a 
time-dependent, but spatially homogeneous, anisotropic 
stress can be evolved in an anisotropic universe to first 
order. If we attribute this anisotropic stress to anisotropy 
in the dark energy which only comes to dominate a late 
times, then we can deduce strong limits on the anisotropy 
of the equation of state of the dark energy. Quantifying 
this in terms of Aw we find that \Aw\ < 2.1 x 10 -4 for 
fl m < 0.5 

Previous work on anisotropic universes: We note that 
Campanelli et al have taken a similar approach in or- 
der to constrain the ellipticity of the surface of last 
scattering. In their calculation they introduce, ad hoc, 
the anisotropic stress associated with an elliptical sur- 
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FIG. 4: Quadrupole likelihood function for a = (solid line), 
0.1 (dotted line), 0.5 (short-dashed line) and 1.0 (long-dashed 
line). The value of C\ has been arbitrarily set to unity. 



face of last scattering, <7y = — ej ec <5j3<5j3 and constrain 
e rcc < 1CP 2 . Our calculation is similar in approach, but 
the anisotropic stress in our model is dynamical and cal- 
culated from the Einstein and conservation equations. 

Our calculations of the sample variance are somewhat 
at odds with those presented in refs. 0, [23j]. There, 
they found a non-zero anisotropic stress caused by an el- 
liptical Universe could decrease the value of the observed 
quadrupole. Their approach was different, however, in 
that for each realization of the stochastic component they 
minimized the value of Ci by changing the orientation of 
the ellipsoid. This is equivalent to assuming that the 
two components are correlated. In our case, we assume 
the two components are uncorrelated, which leads to the 
increased expectation of C^. It could be that there is 
correlation between the two components in our model 
generated at second order and this is presently being in- 
vestigated. It is, however, likely to be a sub-dominant 
effect. 
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